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ABSTRACT 

CD Context. We introduced the technique of model atmosphere calculation with polarized radiative transfer and magnetic line blanketing. However, 
the calculation of model atmospheres with realistic magnetic field configurations (field strength and angle defined relative to the atmosphere 
1 plane) has not been previously attempted. 

, Aims. In the last paper of this series we study the effects of the magnetic field, varying its strength and orientation, on the model atmosphere 
■ structure, the energy distribution, photometric colors and the hydrogen Balmer line profiles. We compare with the previous results for an 
isotropic case in order to understand whether there is a clear relation between the value of the magnetic field angle and model changes, and to 
, study how important the additional orientational information is. Also, we examine the probable explanation of the visual flux depressions of 
["■■■«. ■ the magnetic chemically peculiar stars in the context of this work. 

Methods. We calculated one more grid of the model atmospheres of magnetic A and B stars for different effective temperatures (r c(f = 8000 K, 
*0 1 1 000 K, 15 000 K), magnetic field strengths (B = 0, 5, 10, 40 kG) and various angles of the magnetic field (Q. = - 90°) with respect to the 
CD ' atmosphere plane. We used the LLModels code which implements a direct method for line opacity calculation, anomalous Zeeman splitting of 
pH , spectral lines, and polarized radiation transfer. 

Q . Results. We have not found significant changes in model atmosphere structure, photometric and spectroscopic observables or profiles of 
hydrogen Balmer lines as we vary the magnetic field inclination angle SI. The strength of the magnetic field plays the main role in magnetic 
line blanketing. We show that the magnetic field has a clear relation to the visual flux depressions of the magnetic CP stars. 
Conclusions. We can use the approach introduced in the previous paper of this series, which neglects anisotropy effects, to calculate model 
atmospheres with magnetic line blanketing. This technique seems to be reliable, at least for homogeneous atmospheres with scaled solar 
abundances. 
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1. Introduction The magnetic field influences the formation of spectral 

lines due to the anomalous Zeeman effect. Spectral lines split 
In the previous papers of this series (Khan et al. 2004; into a number of components, increasing the line absorption 
Kochukhov et al. 2005, Paper I; Khan & Shulyak 2006, and the line equivalent width, which is known as the magnetic 
Paper II) we studied the effects of magnetic line blanketing on intensification effect. Hundreds of thousands of lines modify 
model atmospheres of magnetic chemically peculiar (CP) stars, the line opacity, and then change the model atmosphere struc- 
We showed that models with magnetic line blanketing can pro- ture (due to the backwarming effect), as well as modifying the 
duce the characteristic observable features which are typical outgoing radiation. The specific line profile in the presence of 
of the magnetic CP stars such as flux re distribution from the the magnetic field depends on the anomalous Zeeman splitting 
UV region to the visual jLeckronelll973l) and flux depression pattern, the magnetic field strength and the inclination angle 
around 5200 A dKodairalll969h. In our earlier work we made between the magnetic field vector and the light propagation di- 
some assumptions that simplified the modelling technique but rection (often named the line of sight). The strength of the mag- 
still allowed us to carry out a systematic analysis of the role of netic field and the splitting pattern determine the wavelength 
the magnetic field in the line blanketing and its influence on the shift of the split components, while the shape of component 
observables. profiles is defined by the inclination angle and the splitting pat- 

tern via solution of the radiative transfer equation for polarized 
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In a magnetic plane-parallel stellar atmosphere, the light 
propagating in different directions has different angles to the 
magnetic field vector, and thus carries different amounts of en- 
ergy. However, one can expect that the overall effect of the 
enormous number of split spectral lines on line blanketing de- 
pends mainly on number of these components and the mag- 
netic field magnitude. The point is that for different patterns 
the maximum of the ma gnetic intensification occurs for dif- 
ferent inclination angles (IStift & Leonell2003l Sect. 3.1), and 
particular angle dependent variations of split line profiles may 
appear to be insignificant on the line blanketing after averag- 
ing through the solid angle integration procedure for millions 
of components. Consequently, the approach which neglects the 
dependance of the inclination angle of the magnetic field over 
different pencils of radiation and adopts the same angle for all 
of them looks adequate. Such an approach was used in Paper I 
and Paper II. 

Looking back, we can say that Paper I revealed the general 
effects of the enhanced line blanketing due to Zeeman split- 
ting of spectral lines on the backwarming effect, and resulting 
effects in producing signatures on the outgoing radiation. In 
Paper II we introduced polarized radiation transfer in the mag- 
netic model atmospheres and developed the technique of such 
model atmosphere calculations. We showed that neglecting of 
effects of polarized radiative transfer overestimates the influ- 
ence of the Zeeman splitting on the magnetic opacities, and 
thus can not be used as an accurate approximation. 

One aspect that remains to be included for absolutely cor- 
rect magnetic line blanketing modelling is to take into account 
variation of the inclination angle over different directions in the 
stellar atmosphere. This will allow us to conduct a comprehen- 
sive analysis of the approximation adopted in Paper II to find 
probable systematic errors such as over- or underestimation of 
the line blanketing effect. 

In this final paper of the series we compute a grid of LTE 
ID model atmospheres of A and B stars with the same fun- 
damental parameters used in previous works but replacing the 
simplified approach in which the magnetic field is assigned 
with respect to the light propagation direction. Here we specify 
exactly the strength and the angle of the magnetic field with re- 
spect to the atmosphere plane, and vary them across the funda- 
mental modelling parameters space. We then study the result- 
ing effects on the magnetic line blanketing and the consequent 
changes in the model atmosphere structure and the emergent 
flux. 

The paper is organized as follows. In Sect.[2]we describe 
new aspects of the magnetic model atmosphere calculation. 
Sect.[5]presents numerical results and extensive analysis of the 
flux distribution and Sect.[4]summarizes results of this work. 

2. Calculation of magnetic model atmospheres 

The technique of model atmospheres calculation with polar- 
ized radiative transfer and magnetic line blanketing due to the 
anomalous Zeeman splitting is presented in Paper II in detail. 
In this work we fully implement this technique, dropping an 
assumption about isotropic properties of the stellar atmosphere 
penetrated by the magnetic field, which was adopted in Paper I 
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Fig. 1. Reference frame of plane-parallel atmosphere. The z- 
axis is parallel to the unit vector n and the xy plane coin- 
cides with the plane of an atmosphere. The z'-axis (propaga- 
tion direction of light) has an inclination 6 to the unit vector 
n and azimuth <p measured counterclockwise from the x-axis: 
< < 7r , < ip < 2n. The magnetic vector B is located in 
the xz plane, has an inclination y to the propagation direction 
and an inclination O (0 < £2 < n) to the unit vector. The x' and 
y 1 axes are not shown in figure because their choice does not 
modify the main equations used in this paper. 

and Paper II. Model atmospheres calculated in those papers 
were characterized only by the magnetic field strength, and as- 
sumed the same angle y - n/2 between the pencil of radiation 
and the magnetic field vector, for every light propagation direc- 
tion. The models computed in this paper are explicitly defined 
not only by the strength of the magnetic field but also by its in- 
clination angle O to the vertical (see Paper II, Sect. 2.5). Thus 
the angle y depends on the radiation direction, introducing the 
anisotropy to the modelling, and we have to modify the inte- 
gration scheme over solid angle in order to deal with this. For 
further explanation see Fig. ^ (which is the same as Fig. 2 in 
Paper II). 

In the general case, the integration procedure is performed 
over inclination angle 8 to the vertical (0 < 9 < n) and az- 
imuthal angle p (0 < p < 2n). In a case of an isotropic medium 
(the intensity does not depend on azimuth) the integration is 
performed over fi = cos 6 within the limits < yu < 1 instead 
of -1 < fi < 1, splitting the integral into two for outward and 
inward intensities. 

For model atmospheres with magnetic line blanketing it is 
also possible to split the integral and reduce the coverage of /j. 
to the interval < fi < 1, accompanied by the reducing of the 
interval for the ip angle to limits < tp < n. This follows from 
several properties of the symmetry in the magnetic case. The 
first one is that for inward radiation the rotation of the model 
atmosphere reference frame xyz about the x-axis by n (the in- 
clination angle y of the magnetic field to the pencil of radiation 
converts to the angle n -y) does not change the total intensity 
/ of the beam. The second and third ones arise from the expres- 
sion for cos y (see Eq. (46) in Paper II) 

cosy = cos 9 cos Q + sin 9 sin Q cos p, (1) 

which is illustrated in Fig. |2 One can see the mirror symme- 
try for cos y with respect to the line ip - n and the symmetry 
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Fig. 2. The dependance of cos y on the azimuthal angle <p for 
different inclination angles 0. The angle Q. of the inclination of 
the magnetic field is set to 45°. The figure illustrates the sym- 
metric properties of the cosy values. The solid lines represent 
curves for < < n/2 , the dashed lines for n/2 < 9 < n. 

of rotation around the point (n/2; 0) for < ip < n . In other 
words, the xz plane which contains the magnetic vector B nat- 
urally divides both the top (outward radiation) and the bottom 
(inward radiation) integration domains (hemispheres) into two 
equal integrals. At the same time, similar to the isotropic case, 
the integration over the inclination 6 is also reduced by half. In 
fact, the change of the angle yion-y and general rotation of 
the local propagation reference frame (x'y'z') about the z'-axis 
do not change the flux or radiative equilibrium equations, so 
Eqs. (42) and (44) in Paper II hold directly in this study. 

Consequently, the integration over a solid angle is carried 
out on a 2D (0, <p) grid for < < n/2 and < tp < n that 
corresponds to a quarter of a sphere. The integration is per- 
formed over 4 quadrature points in fi — cos 8, and over 8 points 
in azimuthal angle <p, since the range of tp is twice large as that 
of 0; in total we use 32 points. The numerical integration em- 
ploys Gauss-Legandre quadratures. Our numerical tests show 
that this number of points is absolutely sufficient for the inte- 
gration procedure. Even an integration using 3 by 6 points pro- 
vides an accuracy within 1% . This conclusion coin cides with 
the result of extensive tests bv lAlecian & Stifd J2Q02i Sect. 5.1). 

3. Numerical results 

In order to be consistent with previous studies we have used 
the same fundamental parameters as in Paper II. The stellar 
effective temperature assumed has the values T e ff = 8000 K, 
11000K, 15 000 K, the surface gravity logg = 4.0, metallic - 
ities [M/H] =0.0, +1.0 (for non-solar metallicities, the He 
abundance was decreased to normalize the sum of all abun- 
dances), magnetic field strengths 0, 5, 10, and 40 kG, with in- 
clination Q to the atmosphere plane Q = 0°, 45° and 90° (for 
two subsets of models the angle Q ranges from 0° to 90° in 15° 
steps). The parameter values that we have assumed approxi- 
mately cover the stellar parameter space occupied by the mag- 
netic Ap and Bp stars. 



We used the logarithmic continuum optical depth scale 
logT5ooo running from +2 to -6.875 and subdivided into 72 
equally-spaced layers as a scale for the model atmospheres 
calculation. Convection and microturbulent velocity were ne- 
glected because of the presence of a global magnetic field and 
direct inclusion of the magnetic intensification. The wavelength 
ranges for the spectrum synthesis are as specified in Paper II 
(Sect. 3). 

For model atmosphere calculatio n we used a li ne list ex- 
tracted from the VALD database (Kuo ka et alJl9 99) including 
lines originating from predicted energy levels. The procedure 
of the preparation of the line list and dealing with Lande factors 
is described in Paper II (Sect. 3.1). Here we recall that Zeeman 
splitting is taken into account for all spectral lines except the 
hydrogen lines. 

We then used this line list for line preselection in the 
LLModels code, using the selection threshold 1 % (see Paper II, 
Sect. 3.1). For each set of T e g, [M/H], B values, we have used a 
corresponding converged model calculated in Paper II as an ini- 
tial guess of the model atmosphere structure and for the spectral 
line preselection procedure. 

In the following sections we study the model temperature 
structure, spectral energy distribution, photometric colors and 
profiles of hydrogen Balmer lines. Our main interest is to com- 
pare them to results obtained in the previous paper, i.e. to model 
atmospheres calculated assuming that the orientation of the 
magnetic field is characterized by the same angle y = n/2 for 
every light propagation direction. In those models we neglect 
the anisotropy induced by the magnetic field, because the mag- 
netic field is specified with respect to the local directions of 
the pencil of radiation, not to the atmosphere plane. For fu- 
ture reference we will call these model atmospheres "reference 
magnetically isotropic models". 

To perform a self-consistent study and eliminate possi- 
ble systematic errors, we recalculated reference magnetically 
isotropic models because a new flux integration scheme was 
introduced and a new version of the LLModels code was is- 
sued. We also reduced the number of model atmospheres to be 
computed in comparison to our usual grid: we did not consider 
models with intermediate metallicity value [M/H] = +0.5 and 
magnetic field strengths B — 1 kG and 20 kG. This allowed us 
to save computational time without losing a view of the overall 
picture. 

We note that the new models are pretty time consuming; 
one iteration takes about 7 hour on a modern PC or laptop, in 
comparison to about 2 hours and 50 minutes for models calcu- 
lated in Paper II and Paper I respectively, while one iteration 
for non-magnetic models takes about 5 minutes. To get con- 
verged the models presented in this paper requires 6-12 iter- 
ations. In Paper I and Paper II the number of iterations was 
7-12 and 9-20 respectively, depending on modelling parame- 
ters. We recall that in this series we used simplified or previ- 
ously calculated model atmospheres (with the same r e ff, logg 
and [M/H] or even with the same magnetic field strength) as 
an initial guess of the model atmosphere structure. This pro- 
vides us a very good starting point to iterate a new model with 
the more sophisticated approach. Usually, for new models, five 
iterations are enough to achieve convergence throughout the 
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whole atmosphere except surface layers where the rest of the 
iterations are required to correct the temperature of these lay- 
ers by just 2-10 K. Poor convergence in surface layers is a 
known issue since ATLAS9, and controlled by the condition 
of radiative equilibrium (Paper II, Sect. 2.3). Essentially, di- 
rect consideration of the anomalous Zeeman splitting increases 
the number of spectral lines (jt- and cr-components) involved 
in calculations by a factor of about 17. Moreover, the integra- 
tion is carried out on a 32 point 2D grid instead of the usual 
3 yu-point ID grid, and the solution of the polarized radiative 
transfer equations takes more time than the unpolarized one. 
The computation time does not grow linearly with the number 
of integration points since a lot of basic physics involved into 
calculations (like the very time consuming sums of Voigt and 
Faraday- Voigt profiles, see Eqs. (17) and (18) in Paper II) needs 
to be computed only once for all integration points at a given 
wavelength. 

With the same goal of saving computational time, we lim- 
ited ourselves to three values of the inclination angle Q. = 0°, 
45° and 90° for each model of the grid. However, before doing 
this we checked for two different r e jf and B values that the vari- 
ation of the Q angle on a fine angle grid produces smooth and 
consecutive changes on the model atmosphere structure, and 
Q. = 0° and 90° correspond to boundaries of these changes. 

An important point about model atmospheres with mag- 
netic line blanketing that we must discuss to avoid confusion, 
is determination of the stellar flux. The flux that comes from 
a star to the observer from the sum of all visible stellar sur- 
face unit areas, we call the stellar flux. This is exactly the flux 
that observations provide us as energy distribution or photo- 
metric colors. In order to compare theoretical predictions with 
observed values we need to compute the stellar flux. 

In a isotropic medium (absorption coefficient does not de- 
pend on direction) such as a non-magnetic plane-parallel atmo- 
sphere, the intensity of the radiation depends only on its angle 
fj. = cos to the vertical. At the same time, for uniform spheri- 
cal stars, the variation of the intensity across the visible stellar 
disk has the same dependence as the local angular variation, 
and we can cover the stellar surface with the identical mod- 
els. Thus, the calculation of the surface flux per unit area in 
a plane-parallel atmosphere automatically produces the stellar 
flux, which can be compared with observations. 

In an anisotropic medium, such as a medium penetrated by 
a magnetic field, the extinction depends on the light propaga- 
tion direction, thus the flux from a point on the atmosphere sur- 
face does not necessarily equal the stellar flux. In the general 
case of a non-uniform stellar surface, the star has to be cov- 
ered by a grid of different model atmospheres, i.e. as applied to 
our case models need to be calculated with individual values of 
the strength and inclination of the magnetic field. Then, calcu- 
lating the intensity of radiation coming from each local model 
to the observer, and performing the integration over the visible 
spherical surface in the observer reference frame, we obtain the 
stellar flux. 

Since the real configuration of a global stellar magnetic 
field is pretty complicated, and its strength can vary substan- 
tially over the stellar surface, it was difficult to restrict our study 
to certain magnetic field configurations. In this paper our fo- 



cus is to study the magnetic line blanketing effect as a func- 
tion of the strength and inclination of the magnetic field in a 
plane-parallel atmosphere. The modelling of complicated con- 
figurations of the magnetic field naturally would bring us to 
modelling of particular objects, accompanying by much com- 
putational expense, and would not illuminate clearly the basic 
effects of magnetic line blanketing on the outgoing stellar radi- 
ation. 

Thus, we decided to consider the surface flux produced by 
the atmosphere model as a stellar flux. This approach corre- 
sponds to the case in which the magnetic field has the same 
strength over the stellar surface, and the angle between the ob- 
server's line of sight and the local magnetic field vector varies 
in the same way over the stellar surface as the angle y between 
a local pencil of radiation and the vector of the magnetic field 
in a plane-parallel atmosphere. For example, for Q = 0°, i.e. 
for a magnetic field that is perpendicular to the atmosphere 
plane, the surface flux produced by the model atmosphere cor- 
responds to the radial magnetic field configuration in which the 
magnetic field vector is perpendicular to the stellar surface at 
each point. As another example, for Q = 90°, the case in which 
the magnetic field vector lies in the atmosphere plane, the sur- 
face flux produced by the model atmosphere corresponds to a 
case in which the magnetic field vector is a tangent to the stel- 
lar surface at each point. Obviously, that two these cases are 
quite different and should give us an idea what kind of changes 
we may expect in energy distribution. For other Q. angle values 
between 0° and 90° it is rather difficult to draw a mental pic- 
ture of the stellar magnetic field configuration. Some ideas of 
possible fields can be obtained from Fig. [2 

Note that considerations about the stellar flux outlined 
above have nothing to do with the model atmosphere calcu- 
lation technique or model atmosphere structure, and concern 
only interpretation of the spectroscopic observables. 

3.1. Model structure 

In Fig.[3]we show the difference between the temperature struc- 
ture of magnetic angle dependent, and reference magnetically 
isotropic models, as a function of the continuum optical depth 
log 

T 5000 for metallicity [M/ H] — + 1 .0. 

We note that it is difficult to analyse the temperature 
changes using the log T5000 scale, in view of the backwarming 
effect, since the monochromatic continuum depth scale T5000 
does not contain frequency averaged mean absorption coeffi- 
cient (unlike the Rosseland T ross scale). But having model at- 
mospheres computed in log r ross scale (Paper I, Fig. 1) and re- 
plotting the temperature distributions for them on the log T5000 
scale, we can say that the temperature trends are pretty similar 
in both scales from the top of the atmosphere down approxi- 
mately to log Tross = log T5000 = -2. This comparison will help 
us shortly. 

Though, in general, effects on the model structure in 
this work are small, two interesting features can be pointed 
out. The first one is that some essential properties in the 
model atmosphere structure appear only for depths deeper that 
log T5000 = — 1, while the rest of the atmosphere upward from 
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Fig. 3. The temperature differences between magnetic angle dependent and reference magnetically isotropic model atmospheres 
for the effective temperatures r eff = 8000 K, 11 000 K, 15 000 K, metallicity [M/H] = +1.0, strengths of the magnetic field B = 5, 
10, 40 kG, and angles of the magnetic field 0°, 45° and 90°. Two frames for T eff = 8000 K, B = lOkG and T eft - = 15 000 K, 
B = 40 kG represent the same dependence for a fine set of angles ranging from 0° to 90° in 15° steps. The depth scale log T5000 is 
a logarithmic continuum monochromatic optical depth scale calculated at 5000 A. 



this depth is almost unchanged. Taking into account the fact 
about the similarity between changes of the temperature dis- 
tribution in log T5000 and log t toss scales for the optically thin 
upper region, we suppose that the variation of the inclination 
angle O does not have a substantial influence on the line blan- 
keting (and thus on the backwarming effect), since there are 
not any characteristic distinctions between temperature distri- 
butions for different O angles in upper layers (Fig.[3j. 

The second point is that the model structure in the deeper 
layers demonstrates a clear dependency on the variation of the 
inclination angle Q of the magnetic field. The larger the an- 
gle is, the hotter the atmosphere is. The fine angle grid subsets 
of models for r eff = 8000 K, B = lOkG and T eS = 15 000K, 
B = 40 kG show this orderly behaviour in detail in Fig. [3] Note, 
that the model structures for Q. = 0° and 15° are absolutely in- 
distinguishable. That is an interesting point since increasing of 
the Q angle beyond the 15° value produces clear changes in the 
model atmosphere structure. It looks like this behaviour is con- 
nected with the dependance of the magnetic intensification on 
the field inclination angle. 

The solar metallicity ([M/H] = 0.0) models show the same 
features and temperature behaviour as discussed above but the 
effect in model structure is approximately four times smaller 



than for overabundant ([M/H] = +1.0) models, and is thus 
completely negligible. 

3.2. Energy distribution 

We performed a comparison of the energy distributions of mod- 
els with different Q. angles for each combination of the consid- 
ered modelling parameters T e e, [M/H], and B, relative to the 
reference magnetically isotropic models. We look particularly 
for changes in features which were found to be sensitive to the 
magnetic field in the previous papers: the flux redistribution 
from the ultraviolet region to the visual, and flux depression 
around 5200 A. 

We found that the model atmospheres calculated with dif- 
ferent value of the angle Q produce fluxes that are very close 
each to other and to the reference model flux. The discrepancy 
in the UV region is slightly more recognizable than in the visual 
and can reach as much as only 0.04 mag (on the log/^i scale) 
at some UV wavelength points, but no common trends were 
found. For a strong magnetic field B = 40 kG and overabun- 
dant composition [M/H] = +1.0, the reference magnetically 
isotropic models with effective temperatures r e jf = 8000 K and 
1 1 000 K produce UV flux that is a little bit less deficient than 
angle dependant models with the same fundamental parame- 
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ters, so a small offset appears. The typical value of this offset 
amounts approximately to 0.06 mag for r e ff = 8000 K and to 
0.05 mag for r eff = 11000K. 

3.3. Colors 

In order to perform a complementary quantitative analysis of 
possible flux changes in the visual region we calculated photo- 
metric colors using energy distributions. In this work we con- 
sider the Stromgren-Crawford uvbyUB system, the peculiarity 
index a, and the peculiarity indices V\ - G and Z which are 
linear combinations of the Geneva photometric indices. The 
meanings of these indices, corresponding references and a de- 
scription of a computational procedure can be found in Paper I 
(Sect. 3.3). 

We note that we have not shown any dependency of the 
photometric indices with variation of the angle Q., primarily 
because such variations are rather small. For the same reason a 
graphical representation of color values is not illuminating, and 
we provide a short description only. 

For color analysis we consider two difference values de- 
fined among models with different angle Q. while the rest of 
modelling parameters remain intact. The first one we call de- 
viation; it equals the difference between maximal and minimal 
index values. The second one we call the offset; it equals the 
maximum color change with respect to the reference magneti- 
cally isotropic model. Clearly, the absolute value of the maxi- 
mum offset includes a deviation value. 

We found that the deviation does not exceed 1 mmag for 
peculiar indices {a, Z and V\ - G), and 3 mmag for all indices 
of the uvbyHB Stromgren-Crawford system. 

For the extreme case of a strong magnetic field (B = 40 kG) 
and overabundant chemical composition [M/H] = +1.0 the 
deviations above are increased by 1 mmag for the a and 
mi indices for T e g - 8000 K, and for the c\ index for 

r eff = iioook. 

The color offset for all indices and all calculated models is 
usually not more than +3 mmag, with no noticeable trends, so 
it is very difficult to say that there is any systematic differences 
in comparison to the reference model colors, although some 
exceptions appear. 

For the metallicity indicator index m\ for models with 
r eff = 8000 K and [M/H] = +1.0 the offset grows from 
+5 mmag for B = 5kG up to +20 mmag for B = 40 kG. 
Another exception appears for the peculiar a index for 
r eff = 8000 K, B = 40 kG model: the offset of the a index for 
[M/H] - 0.0 is +5 mmag and for [M/H] = +1.0 is +13 mmag. 
For Z and c\ indices exceptions exist for [M/H] = +1.0, 
B = 40 kG models: the Z index shows the offset -6 mmag for 
T e « = 8000 K, and the offset of the c\ index reaches value 
-7 mmag for r eff = 15 000 K. 

It is clearly seen that almost all of the exceptions occur only 
for the extremely strong magnetic field value 40 kG and over- 
abundant chemical composition [M/H] = +1.0. 



3.4. Hydrogen line profiles 

Since hydrogen lines are often used for determination of the 
stellar fundamental parameters we checked the influence of 
changes in the model atmosphere structure on the hydrogen 
Balmer lines. We calculated profiles of Ha, HB and Hy and 
compared them to the corresponding profiles computed for the 
reference magnetically isotropic models. 

We found that the variation of the angle Q produces dif- 
ferences in hydrogen profiles that are almost negligible with 
respect each to other, as well as with respect to profiles of mag- 
netically isotropic models. The common deviation value for the 
overwhelming majority of models is less than 0.1-0.2% of the 
continuum level and reaches the maximum level of 0.4-0.5% 
for r eff = IIOOOK, [M/H] = +1.0, B = 40 kG model for the 
HB and Hy lines. 

We note that hydrogen line profiles calculated within the 
framework of this work for non-magnetic models (for a cross 
check between papers as our code evolves) differ from those 
calculated in Paper II by less than by 0.1%. 

3.5. Stellar atmospheric parameters 

Since the overall modification of the photometric indices of the 
uvbyHB system found in this work ("Sect. 13.31 is quite small we 
do not consider the influence of the angle Q. variation on the 
photometric determination of the stellar atmospheric parame- 
ters. 

However, the insensitivity of the photometric indices to the 
inclination of the magnetic field allows us to reconsider results 
found in Paper II with more confidence. In the previous papers 
we found that the magnetic opacity does not produce system- 
atic errors in the photometric determination of the atmosphere 
parameters of CP stars. The T e R and logg values of the mag- 
netic models derived from the photometric calibration for nor- 
mal stars differ from those for non-magnetic models within the 
usual error bars. It appears that the flux redistribution from the 
UV to the visual region and flux depression features do not af- 
fect essentially the general photometric dependencies used for 
stellar atmosphere parameters determination. We suppose that 
this happens mainly because the slope of the Paschen contin- 
uum does not suffer much from the energy redistribution phe- 
nomenon. The energy excess in the whole visual region ap- 
pears to be evenly shifted along the vertical axis only (Paper II, 
Fig. 3). 

3.6. Visual flux depression 

Having the results obtained in this paper, proving that the incli- 
nation of the magnetic field is not a crucial aspect for construc- 
tion of model atmospheres with magnetic line blanketing, we 
can now explore the fluxes of Paper II with more confidence. 

The spectra of magnetic CP stars show three broad band 
visual flux dep ressions around 4200 A, 5200 A and 6300 A 
(Kodaira 1969). Some authors indicate slightly different central 
wavelengths of these features (e.g. 4100 A and 5300 A) which 
is not surprising because the depression around 4200 A as well 
as around 6300 A spans several hundred A, while the depres- 
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Fig. 4. Synthetic fluxes in the 4200 A region normalized at 5556 A for the effective temperatures r eff = 8000 K, 1 1 000 K, 
15 000 K, metallicities [M/H] - 0.0, +0.5, +1.0, and the magnetic field strengths B = 0, 10, 40 kG. Original stellar fluxes are 
convolved with a Gaussian profile with FWHM = 15 A. 



sion around 5200 A is almost for 1000 A wide (see lAde lman 
ll979tlAde lman & P voeJ 19931 and the whole series). 

The 5200 A feature is detected in most magnetic CP stars, 
the 4200 A feature does not appear for all of them, and only 
a small number of CP stars show the 6300 A flux feature. The 
presence of all these depressions tends to be observed only for 
magnetic CP stars. 

To study these features, the classic photometry systems 
such as Geneva or Stromgren (Hauck 1974: | Masana et all 
1998) and the specialized a system (Maitzen 1976), as well as 
spectrophotometric scans of the energy distribution, have been 
used. 

The origin of the flux depressions are most likely con- 
nected with the line blanketing resulted from an abundance 
peculiarity and magnetic field effects. Recent investigations 
show ed that model atmos pheres with individual abundance pat- 
terns (Kupka et al. 2004), or even individual abundances com- 
bined with stratified abundances dShulvak et all2004l) . can ex- 
plain observed energy distribution of Ap stars while mod- 
els with scaled and h omogeneous abundances fail to do this 
lAdelman et al 1199.4. 

Thus, an individual chemical composition and a non- 
uniform distributi on of chemical elements with depth that 
is observed JWade et alJ l200ll) and theoretically explained 



jLeBlanc & MoninTl2004 ) should without doubt be taken into 
consideration for the atmosphere modelling of individual stars. 



Since this series is focused on effects of the magnetic field 
on the line blanketing, and does not study individual objects, 
we did not consider abundances other than scaled to the solar 
ones. This allowed us to build an uniform grid of model atmo- 
spheres to find the main dependencies between the magnetic 
line blanketing and changes in the energy distribution. 

A part of the energy distribution analysis in this series was 
an investigation of the photometric indices. Here we want to 
call attention to the nt\ index. The color-index my of the uvbyHB 
system is often used as a metallicity indicator of chemically 
peculiar stars. Although it was not introduced as a special pa- 
rameter to separate peculiar stars from normal ones, its value is 
determined by three narrow band filters vby which are affected 
by the 4200 A and 5200 A features. As follows from Fig. 4 in 
Paper I, the m\ index has a clear dependance on the magnetic 
field strength for low effective temperature. The results pre- 
sented here in Sect. !3.3l also suggest that m.\ is sensitive to the 
magnetic field. Thus, we decided to expand our energy distri- 
bution study and look for probable depressions around 4200 A 
and 6300 A (the flux feature around 5200 A was well studied in 
Paper I and Paper II). 

Fig. |4] and Fig. [5] show the spectral energy distributions 
for the 4200 A and 6300 A regions, based on the data of 
Paper II. The synthetic fluxes are convolved with a Gaussian 
and normalized at 5556 A (m,\ - \ogF^). This normalization is 
a usual procedure for observed energy distributions, which are 
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Fig. 5. Synthetic fluxes in the 6300 A region normalized at 5556 A for the effective temperatures T e ff = 8000 K, 1 1 000 K, 
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convolved with a Gaussian profile with FWHM = 15 A. 



often normalize d at 5000 A (|Adelmanet all 1 19951) or 5556 A 
llAdelmanlll975l) . The latter wavelength corresponds approxi- 
mately to the mean wavelength 5500 A of the Geneva V pho- 
tometric band. The normalization of the flux at both of these 
wavelengths can be undesirably affected by the nearby 5200 A 
depression which decreases the true continuum level. However 
our purpose was to plot fluxes in a manner suitable for com- 
parison, and we found that normalization at 5556 A provides a 
better visualization without addi t ional fl ux shifts along vertical 
axis. We note that Kup ka et alJ |2004) showed no significant 
differences between these two kinds of the normalization. 

We found that the depression features in spectral regions 
around 4200 A and 6300 A depend on the effective temper- 
ature, metallicity and magnetic field strength. The usual de- 
pendance on the stellar effective temperature, namely that for 
low r e ff the line blanketing produces substantial changes while 
for higher temperatures this effect becomes rather small, takes 
place again. The sensitivity of depressions to the magnetic field 
magnitude depends on metallicity and grows with increasing 
metallicity. We note that the overabundant chemical composi- 
tion itself for non-magnetic models has a positive effect on the 
depression magnitude. 

The maximum deficiency produced by a magnetic field 
B = lOkG compared to the non-magnetic energy distribu- 
tion for the 4200 A region amounts to 20^40 mmag for 
r eff = 8000K and 5-10 mmag for r eff = 11000K; for the 



6300 A region these values are 10-15 mmag for T e ff = 8000 K 
and about 5 mmag for T e ff = 11 000 K. For the higher effective 
temperature of = 15 000 K, differences between magnetic 
and non-magnetic energy distributions almost disappear. 

Since there are no well-established photometric systems 
like the a system to measure flux depressions around 4200 A 
and 6300 A it is difficult to make a definitive conclusion about 
the relation to the observed values. But it is clear that the mag- 
netic field can be an important factor in these spectral regions 
to increase the magnitude of the depression. Also, it should 
be remembered that a specific chemical composition, and es- 
pecially stratified abundances, can result in more pronounced 
changes in the visual energy distributions than a magnetic field 
of moderate values, but that the magnetic field itself can eas- 
ily intensify weak features and make them visible in spectra 
of magnetic CP stars. A comprehensive analysis of model at- 
mospheres with individualized or stratified abundances and the 
resulting fluxes has never been done. While it is obvious now 
that the study of certain CP objects requires a speci fic model 
atmosphere (Kupk a et aTll2004l IShulvak et al1 E004) we want 
to emphasize the necessity to take into account the magnetic 
line blanketing for a detailed study, especially since we showed 
that a non-zero microturbulent velocity is not able to reproduce 
the effects of the magnetic intensification properly (Paper I, 
Sect. 3.6). 
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4. Conclusions and discussion 

This paper ends the series of works concerning the magnetic 
line blanketing in stellar atmospheres and its relation to the ob- 
servable quantities of magnetic CP stars. Using the technique 
introduced in Paper II, we calculated a grid of model atmo- 
spheres of magnetic A and B stars for different strengths of 
the magnetic field and different angles of its inclination to the 
atmosphere plane. We have used this grid to analyze the be- 
haviour of the model atmosphere structure and the observables 
for variety of fundamental modelling parameters, mainly to 
look for any discrepancy between models with different mag- 
netic field inclination angles, and to compare these models with 
previous results (Paper II). 

We have not found any significant changes in model at- 
mosphere structure, photometric indices, the energy distribu- 
tion and profiles of hydrogen Balmer lines that depend on the 
magnetic field inclination angle Q. Moreover we found that 
model atmospheres for different inclination Q angles are al- 
most the same as models calculated neglecting the anisotropy 
effects (Paper II), when the magnetic field direction is assigned 
with respect to the pencil of radiation, not to the atmosphere 
plane. Only for extreme values of the magnetic field, such as 
B = 40 kG, do the results indicate that magnetically isotropic 
models tend to produce less deficient flux redistribution and 
5200 A flux depression than magnetic angle dependent mod- 
els, while there are no significant differences among these. 

It appears that the magnetic field strength is the main fac- 
tor to affect the line blanketing. We conclude that it is clearly 
possible to implement the isotropic approach (Paper II) to cal- 
culate model atmospheres with magnetic line blanketing. This 
approximation looks reliable, at least for the homogeneous at- 
mospheres with scaled to the solar abundances considered in 
our studies, and can decrease the computation time dramati- 
cally. For atmospheres with stratified or individualized abun- 
dances this conclusion has still to be checked. At the same time 
in this paper we found that variation of the scaled abundances 
clearly has a more pronounced effect on model atmospheres 
than variation of the Q angle. 

For a particular magnetic field configuration (like a combi- 
nation of dipole, quadrupole and octupole) the strength of the 
magnetic field can vary substantially over the visible stellar sur- 
face and individual model atmospheres have to be calculated to 
cover it. However, taking into account a very weak dependence 
of the model atmosphere structure on the magnetic field incli- 
nation to the atmosphere plane found in this study the compu- 
tational expenses can be decreased. The point is that there are 
always unit domains with the same (or almost the same) mag- 
netic field strength on the surface integration grid. Thus, the 
same model atmosphere (computed only once) can be used to 
calculate the outward radiation intensities for a number of sur- 
face units which have very close values of the magnetic field 
strength but different orientations to the observer line of sight. 
Then, integrating individual intensities over all unit domains 
we obtain the stellar flux corresponding to a particular mag- 
netic field configuration. 

Also we note that in general, the magnetic fields of Ap 
stars are considered to be dipoles to first order. The field in- 



tensity for a dipolar field is twice as strong at the poles as 
at the equator. Essentially all studies using longitudinal field 
and surface field measurements have found that the fields de- 
viate from a dipole - that one pole is stronger than the other 
(that a quadrupolar component is required) and often that the 
pole-to-equator variation of the field strength is not as strong 
as that of a dipole (an octu polar component is required) (e.g. 
lLandstreet & Mafhvsl EoOO). Consequently, a typical range of 
field intensity variation over the surface of Ap stars is probably 
about 30-40%. Thus, we suppose that a model atmosphere cal- 
culated using the isotropic approach and the observable mean 
surface magnetic field strength B s as an input model parame- 
ter seems to be reliable (say, for routine abundance analysis of 
magnetic Ap stars). 

Furthermore, for the isotropic approach with y = jt/2 (i.e. 
the magnetic field is perpendicular to every pencil of radiation) 
and with \ — (the azimuth of the magnetic field in the local 
reference frame x'y'z' , and its choice is free, see Fig. [0, the 
line propagation matrix O (Eq. (7) in Paper II) is 



' 4>i 4>q 
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<t>Q <Pl 











4>i 


>ffQ 


,0 




4>i , 



<D 



Thus the system of four Stokes transfer equations 

»^-=K(I-S), 
dm 



(2) 



(3) 



where / = (/, Q, U, V) + is the Stokes vector, m is column mass, 
S = (5,0,0,0)' is the source vector, and K is the total propa- 
gation matrix (see Eqs. (2-8) in Paper II), decouples into two 
independent sets of two equations each one. The second set 
for U and V does not contain emission terms and is out of our 
interest because boundary conditions for inward and outward 
radiation assume initially unpolarized light (see Eqs. (38-39) in 
Paper II). The first set for / and Q is 



dl 

P 'dm = k i( I ~ S "> + k qQ- 
'"dm = k qV ~ + K 'Q> 



(4) 



where kj and kq are the elements of the total propagation matrix 
K which correspond to 4>i and <pQ of the line propagation matrix 
O (see Eq. (3) in Paper II). 

Consequently, the transfer equation for polarized radiation 
reduces to the solution of two standard unpolarized radiative 
transfer equations with different substitutions for the intensity 
and absorption coefficient 



d(I±Q) 
dm 



(k i ±kq)(I+Q-S). 



(5) 



The final Stokes parameters / and Q are defined as a half-sum 
and half-difference of the solutions, respectively. In that way, 
for the magnetically isotropic models with y - n/2, we can 
save even more computational time. 

We also supplemented the investigation of the most promi- 
nent feature of the CP stars in the visual region with the flux 
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depressions centered at 4200 A and 6300 A. Our numerical re- Russell Diagram, ed. G. Mathys, S. K. Solanki, & D. T. 
suits show that the magnitude of these features clearly depends Wickramasinghe, ASP Conf. Ser., 248, 373 
on the magnetic field strength, effective temperature and metal- 
licity. The depressions grow with increasing magnetic field 
strength and metallicity. However, for higher T e g the magni- 
tude and the sensitivity to the magnetic field goes down quickly, 
which is similar to the results for the depression around 5200 A 
(Paper II). 

Taking into account the common relation between vi- 
sual flux depression magnitudes around 4200 A, 5200 A and 
6300 A, and the effective temperature found in this series, we 
suppose that individual chemical composition and inhomoge- 
neous vertical abundance distribution are most probable can- 
didates to explain visual peculiarities of Ap and Bp stars with 
higher effective temperatures. 
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